Asymptotic Statistics of Poincare Recurrences 
in Hamiltonian Systems with Divided Phase Space 
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By different methods we show that for dynamical chaos in the standard map with critical golden 
curve the Poincare recurrences P(t) and correlations C(r) asymptotically decay in time as P oc 
C/t oc 1/t 3 . It is also explained why this asymptotic behavior starts only at very large times. We 
argue that the same exponent p = 3 should be also valid for a general chaos border. 

PACS numbers: 05.45. +b 



During last two decades the local structure of phase 
space in chaotic hamiltonian systems and area-preserving 
maps had been studied in great detail These 
researches allowed to understand the universal scaling 
properties in a vicinity of critical invariant curves where 
coexistence of chaos and integrability goes on smaller 
and smaller scales in the phase space. The most stud- 
ied case is the critical golden curve with the rotation 
number r g = [111...] = (V5 — l)/2 for which the scal- 
ing exponents were found with high precision and it was 
shown that the phase space structure is self-similar and 
universal 0. The most studied map with mixed inte- 
grable and chaotic components is the standard map ||] 
where the golden curve is critical at the chaos parame- 

It is believed that for 
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K > Kg all invariant Kolmogorov-Arnold-Moser (KAM) 
curves are destroyed [|J. 

While the local structure of divided phase space is now 
well understood, the statistical properties of dynamics 
still remain unclear in spite of simplicity of these sys- 
tems. Among the most important statistical characteris- 
tics are the correlation function decay in time C(r) and 
the statistics of Poincare recurrences P(t). The later is 
especially convenient for numerical simulations due to its 
natural property P(r) > and statistical stability. The 
first studies of -P(t) in a separatrix map showed that 
at a large time the recurrences decay as a power law 
P(t) tx 1/t p with the exponent p « 1.5 [Q. Investiga- 
tions of other different maps also indicated approximately 
the same value of p |]|| even if it was remarked that p 
can vary from map to map, and that the decay of P(r) 
can even oscillate with lnr f^-|l0[|. Such a slow decay of 
Poincare recurrences was attributed to the sticking of tra- 
jectory near a critical KAM curve which restricts the re- 
gion of chaos J7|-[T2[. Indeed, when approaching the crit- 
ical curve with the border rotation number n, , the local 
diffusion rate D n goes to zero as D n ~ |r(, — r„|"/ 2 ~ l/q™ 
with a — 5 [|ll| where r n — p n /q n are the convergents 
for rf, determined by continued fraction expansion. The 
theoretical value a = 5 was derived from a resonant the- 
ory of critical invariant curve [ll]Jl2] and was confirmed 



by numerical measurements of local diffusion rate in the 
vicinity of critical golden curve in the standard map p3[ . 
Such a decrease of the diffusion rate near the chaos border 
would give the exponent p = 3 if to assume that every- 
thing is determined by the local properties of principal 
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resonances p n /q n given by convergents of rt, [11 



0- 

However, the value p — 3 is strongly different from the 
numerically found p ~ 1.5. Moreover, the special sim- 
ulations for the standard and separatrix maps with the 
border rotation number r& = r g have given different be- 
havior of P(t) and different p [|7 10 in spite of the fact 
that the local structure of the golden critical curve is 
universal. Different attempts were done to reserve this 
difficulty. In [jl5| the authors argued that a contribution 
from non principal resonances can reduce the exponent 
down to p = 2. Other arguments based on disconnection 
of principal resonance scales were proposed in [jl^| , while 
Murray discussed a possibility that larger times are re- 
quired to see p = 3 decay Jll| . During these years differ- 
ent Hamiltonian systems were studied where the values 
of p w 1 - 2.5 have been found §,0-|o|. 

The analysis of Poincare recurrences is interesting not 
only by itself but also because they are directly related 
to the correlation function of dynamical variables [p[-^5[ : 



C(t) ~ n( T ) ~ tP{t)/ < t >~ l/r p - 



(1) 



Here /u(r) is the probability to remain in a given re- 
gion for time t > r, proportional to its measure, and 
< t > is the average recurrence time. This relation can 
be understood as follows. According to the definition 
P{t) = N T /N where is the total number of recurrences 
and N T is the number of recurrences with time t > r. 
Therefore, for the total time T =< t > N on which 
one investigates the recurrences of a trajectory, we have 
P(t) =< t > N T /T ~< t > /i(r)/r where, due to ergod- 
icity of motion the measure (probability to stay) is pro- 
portional to the ratio of time the trajectory spends in the 
region (T T - tN t ) to the total time T (jjl(t) - T T /T). 
Inside the sticking region the dynamical variables are cor- 
related so that C(t) ~ h(t) Jz|Js| [3~T|] . Since the correla- 
tions are directly related to a diffusion rate (D c ~ J Cdr) 
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the exponent p < 2 can lead to a superdiffusive dynamics 
|j] |l2| . For the standard map such a behavior was indeed 
observed in P, fuj|]2l| . All this show that the asymptotic 
decay of Poincare recurrences is a cornerstone statistical 
problem of two-dimensional maps. 

To understand the asymptotic properties of P(r) we 
used for the first time a new approach based on a direct 
computation of exit times from a vicinity of the critical 
golden curve in the standard map 

y — y — K/(2tt) sin(27rx) , x — x + y modi (2) 

with parameter K = K g . The properties of this curve 
had been studied in great detail |J. In particular, the po- 
sitions of unstable fixed points of resonances p n /cin are 
known with high precision To determine the exit 
time T n from the scale q n we placed 100 orbits in a very 
close vicinity of an unstable fixed point and computed 
the average exit time. For each orbit the exit time is 
determined as a time after which the orbit crosses the 
exit line. The exit line was fixed as y — 1 for the or- 
bits from the side of the main resonance q = 1 or as 
y = 0.5 + asin(27ra;) for the orbits from other side of 
the critical curve with q = 2. In the later case the exit 
line was drown in such a way to cross the two unstable 
points of resonance q — 2 (a = 0.0773...). This allowed 
to take into account the deformation of the q = 2 reso- 
nances. The average exit time t„ from a given scale q n is 
related to the distance of this resonance from the curve 
and is proportional to this distance (measure of chaos) 
Mn = \ r g ~~ r n\ ~ 1/V5Q 2 divided by the local diffusion 
rate D n : r„ ~ ^ n /D n ~ q n . This gives (i ~ C ~ 1/r 2 
and p = 3. The numerical data for dependence of \x n on 
r n is shown in Fig.l. From both sides of the r g curve the 
exit times converge to the asymptotic dependence 



(r 3 /r„) 2 /V5, 



T„ q„ 



T g = 2.11 ■ 10 5 (3) 



This dependence corresponds to the scaling near the 
critical curve Indeed, the local diffusion rate in y 

on the scale q n is D n w AD /q^ , where D = K 2 /8ir 2 
is the quasilinear diffusion rate |22| and A w 0.0066 is 
a numerical constant which is quite small due to a slow 
diffusion inside the separatrix layer [O. As a result the 
sum of transition times between the two scales from r n 
to r„_2 is approximately equal to the total exit time: 
Tji k, J2 n \r n - r n - 2 \ 2 /D n w 1.4 • 10 5 q n - This estimate 
gives the value of r g close to (||) and allows to under- 
stand the physical origin of its large numerical value. It 
is interesting to note that the data of Fig.l show that 
the convergence of r„ to its asymptotic value can be sat- 
isfactory described as \r n /q n T g — 1| oc l/q n . This in- 
dicates a certain similarity between the ratio r n /q n T g 
and the residue R n for periodic orbits which converges 
to its critical value in a similar way [g. Physical rea- 
son of this similarity is the following: R n is related to 
the orbit stability and the larger it is more rapidly the 



orbit escapes from the scale q n . Due to that for odd 
n(q n = 1, 3, 8, ...) the time r„ is smaller than the asymp- 
totic expression (||) (Rn > R cr = 0.250... 0) while for 
even n (q n = 2, 5, 13, ...) it is larger than (|) (R n < R cr 
Q). Due to universality of the critical golden curve struc- 
ture it is natural to expect that the relation ^ and the 
time T g are universal for all area-preserving maps as well 
as R cr (note that q n is the period of unstable periodic 
orbit on the scale n). 

The relation ||) determines the measure of chaotic re- 
gion /i ~ fi n at which a trajectory is stuck for a time 
r ~ T n . Then, according to (Q) the exponent of Poincare 
recurrences is p = 3 and correlations C(t) ~ /i decay as 
inverse square of time. However, this asymptotic decay 
starts in fact only after a very long time r > 10 6 due 
to the large value of T g . This strong delay of asymptotic 
behavior is responsible for the nonuniversal decay ob- 
served for P(t) in (^Jl^l at K — K g . Indeed, for r < r g 
a trajectory doesn't feel the border and [i remains ap- 
proximately constant giving p = 1 that had been seen 
in HJl(]] (see Fig. 2). However, to observe the theoretical 
exponent p = 3 one should go to longer times. To check 
these theoretical expectations we made extensive numer- 
ical simulations of P(t) at K = K g with recurrences on 
the exits lines defined above on the both sides of the crit- 
ical r g curve. The results are presented in Fig. 2 and show 
a clear change in decay of P(r) for r > 10 5 (side = 1) 
and r > 10 7 (side q = 2). To check the relation (|l|) we 
computed P(r) from the data of Fig.l taking the recur- 
rence time r = 2r„ and P(t) — B q < t > q fi n /T where 
q = 1, 2 mark the side of critical curve. With the average 
recurrence time < t >i« 24.5 (or < t >2~ 61) and an 
arbitrary constant B\ = 2.0 (or B>2 = 8.2) the data from 
Fig.l describe the variation of P(t) in the interval of 6 
(4) orders of magnitude. This gives additional support 
for the theoretical exponent p — 3. 

Another check of the relation (0) was done by comput- 
ing the diffusion rate in phase z — z+ (y+y — 2z q )/2 with 
z q = 0(q = 1) and z q = l/2(q = 2). Similar approach 
was used in (21j. The diffusion rate is D c = (Az) 2 /A£ 
and its dependence on time is determined by the decay 
of correlation function of y(t). According to (Q) we have 
D c (t) = D cq G q J tP(t)/ <t > dr = D cq G q J C{t)cLt 
where D cq = \r g - r q \ 2 /3 = 0.049(0.0046) is the quasi- 
linear diffusion rate for q = 1 (2) side and G q ,G q 
are some constants. The correlation C(r) was computed 
from the linear interpolation of data in Fig.l. In addition 
we took that C(2r n ) = /i n //i„i where fj, n i is the value of 
fj, n at the first scale on each side of the critical curve 
with q = 3,5. In this way C(r) remains constant up to 
the first exit time r„i (C(0) = C(r n i) = 1) that corre- 
sponds to the fact that P(t) ~ 1/r for r < r n i ~ r g . 
For p = 3 the rate D c should be finite. The value of D c 
was computed for 100 orbits initially located near un- 
stable fixed point of period q — 1,2. The diffusion rate 
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dependence on r, and comparison with its computation 
from -P(t) and exit times r„ = r/2 via the above in- 
tegral relation, are shown in Fig. 3. Both methods give 
a good agreement with D c (t), especially in the case of 
P(t), confirming |T]). The constants are G\ ~ G2 ~ 2, 
G\ w 0.3, G2 ~ 0.6. According to the above values of 
coefficients -B g the ratio G q /G q should be approximately 
2 times smaller. This may be due to approximate scheme 
used to relate C(r) with ju(t). For r > 10 7 the asymp- 
totic value p — 3 leads to a saturation of D c growth in 
time. Even if the asymptotic diffusion rate is constant 
D c = D c (oo) the distribution function is non-gaussian 
since the higher moments diverge. For smaller r the dif- 
fusion rate D c grows approximately linearly that corre- 
sponds to an intermediate value of p ~ 1. This interme- 
diate slow decay is responsible for the enormously large 
ratio of the asymptotic diffusion rate to its quasilincar 
value: D c (oo)/D cq w 3 • 10 5 (q = f); 10 7 (q = 2). 

The ensemble of data in Figs, f- 3 shows that the 
asymptotic decay of Poincare recurrences and correla- 
tions is determined by the universal structure in the 
vicinity of the critical golden curve, and the contribu- 
tion of boundaries of other internal islands of stability 
is not significant at variance with jf5). The hypothe- 
sis of dynamical disconnection of scales |f(| is also ruled 
out. Our results are in agreement with previous numer- 
ical observations indicated that long recurrences are re- 
lated to orbits being very close to r g curve (|Jl^]. ft 
is interesting to ask the question how the value of the 
exponent p = 3 would be modified for a case of main 
border curve r& = [di, 0,2, Oj, ...] with non-golden con- 
tinued fraction expansion. The numerical data for border 
invariant curves obtained in @] show that the elements 
a, mainly take the values f,2,3,4 and the probability 
to find <Zj > 4 is rather small. For such bounded val- 
ues the general resonant approach developed in 
still show that the diffusion rate near the border scales 
as D n cx l/q n and therefore the exit time will scale as 
T n ~ Qn giving the exponent p = 3. Due to similarity be- 
tween r n /q n and the residue R n discussed above we can 
expect that for a typical the ratio T n /q n will not con- 
verge to a constant as it was for r& = r g but will oscillate 
in a bounded interval similarly to oscillations of R n (see 
e.g. |frof ). Due to above reasons we can expect that even 
in the general case of non-golden border invariant curves 
the average asymptotic universal exponent is p = 3. 

If the exponent p = 3 then it is natural to ask why 
the previous computations of different groups were giv- 
ing p f» 1.5. Our explanation is based on the following 
arguments. First, even for the critical golden curve the 
asymptotic regime starts after very long time which is 
determined by first resonance scales. The first scales are 
not universal that explains why p was varying from map 
to map. If the border curve is non-golden then the ratio 
Tn/Qn should oscillate with n and the asymptotic regime 
will appear even later than for r& = r g . Also on the first 



scales a given map can be often locally approximated by 
the standard map with K w K cr + (r—ri,)df(r)/dr, where 
/ is some smooth function of rotation number |6| . A typi- 
cal example is the separatrix map j6|,p|JTc|] . In this case at 
a given r n the local order parameter is supercritical with 
K — K cr oc \r n — r&l oc l/q n - This scaling is different from 
the asymptotic one with K — K cr oc l/q n oc \r n — r-f, | 1 / 2 
p|, ^0|p6| and can give very long exit time for first scales. 
Indeed, in the standard map with K > K g the transition 
time from y — to y = 1 is proportional to 1/\K — K g \ 3 
HQ that will give exit time t„ ~ j\K - K g \ 3 ~ q^. Ac- 
cording to (0) this will give p — 4/3 that is not far from 
the average p ss 1.5. In addition, when being close to the 
critical curve, as in the standard map with K = K g , one 
should still wait a long time r g to reach the asymptotic 
exponent p = 3. 

In conclusion, we showed that in the case of criti- 
cal golden curve the asymptotic exponent for decay of 
Poincare recurrences is p = 3. This implies that the 
correlation integral converges and the diffusion rate pro- 
duced by such dynamical chaos is finite. However, the 
higher moments of distribution will diverge. We argue 
that the asymptotic exponent should remain the same 
also in the case of a typical border invariant curve. 
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FIG. 1. Dependence of exit time t„ from the scale r„ 
on a distance (chaos measure) from the critical golden curve 
fi n — \r g — r n \ for q n = 3, 8, ...,6765 (black circles) and 
q n — 5, 13, 4181 (open circles). The straight line shows 
asymptotic behavior (3). Error bars are less than the symbol 
size. Logarithms are decimal in Figs. 1-3. 



FIG. 2. Poincare recurrences in the standard map at 
K = K g from the side of resonance q = 2 (upper full curve) 
and q = 1 (lower full curve, shifted down by 3 for clarity). 
Open and full circles show the values of P(r) recalculated 
from the data of Fig.l (see text). The dashed straight lines 
mark the asymptotic decay with theoretical exponent p — 3; 
slope p = 1 is shown by dashed-dotted line. Data for P(r) 
are obtained from 10 orbits of length 10 11 . 



FIG. 3. Dependence of diffusion rate D c on time (full 
curves) compared with its values computed from the Poincare 
recurrences of Fig. 2 (dashes curves) and exit times of Fig.l 
(open and full circles) according to relation (1) (see text). 
Lower curves and circles are for q = 1 side, while the upper 
ones are for q — 2 (shifted up by 4 for clarity). For clarity, 
all circles are shifted up by 0.3 from their optimal positions 
given by coefficients G q (see text) . 
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